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Abstract 

We present a preliminary analysis of the decays K*^^ and K*'^"f using a sample 

of 383 million BB events collected with the BABAR detector at the PEP-II asymmetric energy 
B factory. We measure the branching fractions B{B^ K*^^) = (4.58 it 0.10 ± 0.16) x 10~^ 
and B{B+ K*+'y) = (4.73 ± 0.15 ± 0.17) x 10"^. We measure the direct CP asymmetry to 
be -0.043 < A{B K*-f) < 0.025 and the isospin asymmetry to be -0.021 < Aq- < 0.079, 
where the limits are determined at the 90% confidence interval and include both the statistical and 
systematic uncertainties. 
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1 INTRODUCTION 



In the Standard Model (SM), the decays B — s- K*j [1] proceed dominantly through one-loop 6 — > 57 
electromagnetic penguin transitions. Extensions of the SM predict new high-mass particles that can 
exist in the loop and alter the SM prediction of the branching fractions. The theoretical predictions 
of the decay rates [2~5] for B — > K*'y suffer from large hadronic uncertainties, and previous mea- 
surements of the branching fractions (Table [T]) are more precise than SM estimates. The theoretical 
estimates and experimental measurements of the branching fractions are in reasonable agreement. 

Experimental and theoretical uncertainties are much reduced when considering the CP and 
isospin asymmetries [9], which are defined by: 

^^_ r(B^K_*l)-T{B^K*^) 

T{B ^ K*-f) + T{B ^ K*j) ' ^ ' 



r(^° ^ K*%) - T{B- ^ K*--f) 



(2) 



The K*^ KgTT^ mode is excluded from the determination of the CP asymmetry. Being more 
precise, these quantities allow the SM to be more stringently tested. The SM predictions for the CP 
asymmetry [10] are on the order of 1%, while the isospin asymmetry [5,11] ranges from 2-10%. The 
experimental measurements (Table [I]) are in good agreement with these predictions. However, new 
physics models could alter the SM estimates significantly [11-13], and thus precise measurements 
constrain new physics parameter space. 

This note reports on a measurement of the branching fractions B{B^ K*^j) and B{B~^ — > 
K*~^^), the isospin asymmetry Aq-, and the direct CP asymmetries, A{B^ K*^j) and A{B~^ —^ 





CLEOII [6] 


BABAR [7] 


Belle [8] 






81.9/6-1 


78.0/6-1 


(xlO-5) 


4.55:;j:-±o.34 


3.92 ± 0.20 ± 0.24 


4.01 ±0.21 ±0.17 


B+ K*+'y 


3.76t^|;^ ± 0.28 


3.87 ±0.28 ±0.26 


4.25 ±0.31 ±0.24 


(xlO-5) 






A 


-F0.08± 0.13 ±0.03 


-0.013 ±0.036 ±0.010 


-0.015 ±0.044 ±0.012 


Ao- 


N/A 


+0.050 ± 0.045 ± 0.028 ± 0.024 


±0.012 ± 0.044 ± 0.026 



Table 1: Previous measurements of the branching ratios and asymmetries. The first and second 
errors are statistical and systematic respectively. The last error on the isospin asymmetry for the 
BABAR measurement refers to the error on the production ratio of charged to neutral B events, 
= r(r(45) ^ B+B-)/T{T{AS) B°B^). 



2 THE BABAR DETECTOR AND DATASET 

We use a data sample containing 383 million BB events, corresponding to an integrated luminos- 
ity of 347 fb"! collected at the T{AS) resonance, taken with the BABAR detector at the PEP-II 
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asymmetric-energy e^e~ collider located at the Stanford Linear Accelerator Center (SLAC). These 
results supersede the previous BABMl measurements [7] . 

The BaBAR detector is described in Ref. [14]. Two components that are especially important 
for this analysis are the Csl Electromagnetic Calorimeter (EMC), used to identify and measure 
photon energies, and the DIRC Cherenkov detector, used to identify charged particles. 

3 ANALYSIS METHOD 

We reconstruct K*^'y using the modes K*^ ivT+vr" and K*^ Kstt^, and i?+ K*~^^ 

using the decay modes K*~^ — > iiT+Tr" and K*~^ — > KsTr~^- A high energy photon is combined with 
each vector meson. 

The dominant source of background is continuum events {e'^e~ qq, with q = u,d, s, c) that 
contain a high-energy photon from vr'' or rj decay. The remaining background consists primarily of 
initial-state radiation (ISR) processes, and higher-multiplicity 6 — > 57 decays, where one or more 
particles has not been reconstructed. In addition, the decays of B ^ K*j can enter the signal 
selection by mis-reconstructing a similar mode. For example, the decay B K*j{K*'^ — K+7r°) 
provides background for the mode B K*^{K*^ K~^7t~) by not correctly reconstructing the vr 
meson. For each signal decay mode, selection requirements described below have been optimized 
for maximum statistical sensitivity with an assumed signal branching fraction of 4.0 x 10^^ [7]. 

Photon candidates are identified as localized energy deposits in the EMC that are not associated 
with any charged track. The primary photon candidate is required to have a center-of-mass (CM) 
energy between 1.5 and 3.5 GeV, to be well-isolated and have a shower shape consistent with an 
individual photon [15]. In order to veto photons from 7r° and 7] decays, we form photon pairs 
composed of the signal photon candidate and all other photon candidates in the event. We then 
reject primary photon candidates consistent with coming from a tt'^ or 77 decay based on a likelihood 
ratio that uses the energy of the partner photon, and the invariant mass of the pair. 

The charged tracks must be well-reconstructed in the drift chamber, and are required to be 
consistent with coming from the e"^ e~ interaction region. They are identified as if or tt mesons by 
the Cherenkov angle with respect to track direction, as well as by energy loss of the track (dE/dx). 
The Ks candidates are reconstructed from two oppositely charged tracks that come from a common 
vertex. We require the invariant mass of the pair to be 0.49 < m^+^- < 0.52 GeV/c^(0.48 < 
m^+^- < 0.52 GeV/c^) and have a Ks flight length significance requirement of 9.3(10) for the 
i^*o ^ KsTr'^{K*+ KsTT^) mode. 

We form vr'' candidates by combining two photons (excluding the primary photon candidate) in 
the event, each of which has an energy greater than 30 MeV in the laboratory frame. We require the 
invariant mass of the pair to be 0.112 < m-y-y < 0.15 GeV/c^ and 0.112 < rrijj < 0.15 GeV/c^ for 
the K*^ — > KgTT^ and K*'^ —>■ K'^tt^ modes respectively. In order to refine the 7r° three momentum 
vector, we perform a mass-constrained fit of the two photons. 

We combine the reconstructed K or it mesons to form K* candidates. We require the invariant 
mass of the pair to satisfy 0.78 < mx+Tv- < 1-1 GeV/c^, 0.82 < mj^^^o < 1.0 GeV/c^, 0.79 < 
niK+TrO < 1-0 GeV/c^, and 0.79 < mKg-w+ < 1-0 GeV/c^. The charged track pairs are required to 
originate from a common vertex consistent with the e+e" collision region. 

We combine the K* and high-energy photon candidates to form B candidates. We define in the 
CM frame (the asterisk denotes the CM quantity) = — -Bbeam' where E*^ is the energy of 
the B meson candidate and -Ebeam beam energy. We also define the beam-energy-substituted 
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mass rriEs = y-^beam ~Ps^' where is the momentum of the B candidate. In addition, we 
consider the hehcity angle 9h of the K* , defined as the angle between one of the daughters of the 
K* meson and the B candidate in the K* rest frame. Signal events have AE close to zero with a 
resolution of approximately 50MeV, and an ttt-es distribution centered at the mass of the B meson 
with a resolution of 3 MeV/c^. Since the K* recoils against a photon, it has a cos Oh distribution 
of We only consider candidates in the ranges —0.3 < AE < 0.3 GeV, mss > 5.22 GeV/c^, 

and I cos^j:/! < 0.75. The latter selection is to reject background such as -B ^ K*r] and B K*7r^, 
which are distributed as cos^ ^ in cosOh. To ensure the events are properly reconstructed, we 
apply a selection criterion to the separation (and its uncertainty) along the beam axis between the 
B meson candidate and the rest of the event (ROE) . The ROE is defined as all charged tracks and 
neutral energy deposits in the calorimeter that are not used to reconstruct the B candidate. 

In order to reject continuum background, we combine 13 variables into a neural network 
(NN) [16]. One class of these variables exploits the topological differences between spherical signal 
events and jet-like continuum events by considering information from the B meson candidate and 
the ROE. The other class exploits the difference in particle production mechanisms between B 
meson decays and continuum events. The discriminating variables are described in Ref. [17]. Each 
mode has a separately trained neural network. The output of this network peaks at a value of 
one for signal-like events. We select events with a criterion on this output that is optimized for 
maximum statistical sensitivity. To validate the neural network, we use a B ^ Dn control sample. 

After applying all the selection criteria, we select the best candidate in each event by choosing 
the candidate with the reconstructed K* mass closest to the nominal mass. On average, across all 
four modes, there are approximately 1.1 candidates per event in signal events. 

We perform an unbinned maximum likelihood fit to extract the signal yield, constructing a 
separate fit for each mode. We use three observables (mES) AE, and cos 9h) for each candidate 
event and assume three hypotheses (signal, continuum, and BB) from which the candidate can 
originate. All BB background is included in the BB component. The use of cos 9h suppresses the 
BB background. Since the correlations among the three dimensions are small, we use uncorrelated 
probability distribution functions (PDF) to construct the likelihood function. The correction to 
this method is determined in section HI The likelihood function is: 

/ M \ / ^ 
£ = exp y-Y^Uij ■ ( J] 

where N is the number of events, M is the number of hypotheses, rii represents the yield of a 
particular hypothesis, Vi^xj; oti) is the product of one-dimensional PDFs over the three dimensions, 
Xj = {m-Es, AE,cos9h), and the di represent the fit parameters. 

The signal rriEs PDF for the K*^ K~^tt~ and K*^ Ksi^^ modes is parameterized as 

(3) 

where is the peak position of the distribution, cti^h are the widths to the left and right of the peak, 
and aL^R are a measure of the tails to the left and right of the peak, respectively. We constrain 
(^L = (^R, and fix ai^R to the values obtained from Monte Carlo (MC) simulation [18]. For the 
K*^ KsT^^ and K*^ K~^t^^ modes, the signal mEs distribution is described by a Crystal Ball 
function [19]. The Crystal Ball function has a single tail parameter, a, which we fix to the value 



M 



y^^niVi{xj; 



exp 
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determined from MC. For each mode, the signal IS.E distribution is described by the same function 
in Eq. [31 but with different values for the parameters. However, we allow ai and to float 
independently, but still fix the values of aL,R to MC. For all components, the cosOh distribution is 
modeled by a low order polynomial, which is fixed to the MC values. For the continuum hypothesis, 
the TTT-ES PDF is parameterized by an ARGUS function [20] , with its shape parameter floating in the 
fit. The continuum IS.E shape is modeled by a low order polynomial with its parameters floating in 
the fit. Various functional forms are used to describe the BB background, all parameters of which 
are taken from MC simulation and held fixed. 

The CP asymmetry A parameter is measured in the three "self-tagging" modes: K*^ K~^7r~ , 
K*'^ —^ K~^ir^ and K*~^ Kst^~^ ■ The fit is accomplished by performing a simultaneous fit to the 
two flavor sub-samples {K* and K*) in each mode. All shape parameters are assumed to be flavor 
independent and the A of each component is floated in the fit. 

Figures [1] through H] show the projections of the likelihood fit to data. For each projection, 
signal region cuts (5.27 < ttt-es < 5.29 GeV/c^, —0.2 < IS.E < 0.1 GeV) have been applied, except 
the mES selection is not applied to the rriES distribution and similarly for AE. The asymmetry 
of the signal component of the cos 9h distributions is due to mis-reconstructed signal candidates. 
Table [2] shows the results for the branching fractions and CP asymmetry, where the sign of A is 
defined by Eq. [TJ 



Mode 


6(%) 


Ns 


fi(xlO-^) 


A 


K+TT- 


20.6±0.7 


2394.1±55.6 


4.55 ±0.11 ±0.16 


-0.023 ±0.022 ±0.011 




11.7±0.8 


256.0 ±20.6 


5.01 ±0.40 ±0.37 


N/A 




13.7±0.7 


872.7±37.6 


5.05 ±0.22 ±0.27 


+0.033 ±0.039 ±0.011 




18.8±0.7 


759.1±33.8 


4.56 ±0.20 ±0.17 


-0.006 ±0.041 ±0.011 


B+ K*+-f 
B K*j 






4.58 ±0.10 ±0.16 
4.73 ±0.15 ±0.17 


-0.009 ±0.017 ±0.011 



Table 2: The signal reconstruction efficiency e, the fitted signal yield Ns, branching fraction B, and 
CP asymmetry (A) for each decay mode. The signal efficiencies have been corrected for differences 
between the selection efficiency in data and MC. Errors are statistical and systematic, with the 
exception of e and Ns, which have only systematic and statistical errors respectively. Also shown 
are the combined branching fractions and CP asymmetry. 



4 SYSTEMATIC ERROR STUDIES 

Table[3]lists the sources of systematic uncertainty for all four modes. These are associated with the 
signal reconstruction efficiency, modeling of the BB background, and the choice of fixed parameters 
of the fit PDFs. The "Photon selection" systematic error is a combination of the photon efficiency, 
the isolation criteria, and the shower shape selection. For the Neural Net and the tt^/t] veto, we use 
a S — > Dtt control sample to determine the systematic error. The "Fit Model" systematic error is a 
combination of incorporating uncertainties due to our imperfect knowledge of the normalization and 
shape of the inclusive B Xs'j spectra, and the choice of fixed parameters. We also perform a series 
of experiments in which we select signal events from MC simulation and combine them with events 
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Figure 1: K*^ — > K+vr projection plots of the full fit to data. The daughter of the K* used 
determine the helicity angle is the K meson. 
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Figure 2: K*^ — > Kstt projection plots of the full fit to data. The daughter of the K* used 
determine the helicity angle is the Ks- 



13 




-0.6 -0.4 -0.2 0.2 0.4 0.6 

COSG^ 



Figure 3: K*~^ — > K~^it^ projection plots of the full fit to data. The daughter of the K* used 
determine the helicity angle is the K~^. 
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Figure 4: K*~^ — >■ Kstt~^ projection plots of the full fit to data. The daughter of the K* used to 
determine the helicity angle is the vr meson. 
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from background generated using PDFs from the fit. The bias resulting from correlations among the 
three dimensions, or the PDFs incorrectly modeling the signal distribution can be determined using 
this procedure. The "Signal PDF bias" systematic error results from these series of experiments. 
Associated with all of the systematic uncertainties is a correction factor, which is a ratio between 
the estimated efficiency in data and the corresponding efficiency in MC. The corrections are 0.953, 
0.897, 0.919, and 0.936 for the K*^ K+tt', K*^ Kst^^, K*+ K+7r°, and K*+ Ks7r+ 
modes respectively. We use this factor to correct the MC reconstruction efficiency. 

The systematics of the A measurement were studied in detail in Reference [15]. They were 
found to be due to uncertainties in the hadronic cross section asymmetry and to reconstruction 
asymmetries. Here, we simply adopt the value 1.1%, which is a conservative estimate due to 
reconstruction improvements. 



Table 3: Systematic errors (in %) of the branching fractions. 



Mode 


K*^ ^ K+TT- 


i^*" Kstt" 


K*+ K+tt" 


K*+ KsTT^ 


BB sample size 


1.1 


1.1 


1.1 


1.1 


Tracking efficiency 


1.2 




0.6 


0.8 


Particle identification 


0.6 




0.6 


0.2 


Photon selection 


2.2 


2.2 


2.2 


2.2 


tt" reconstruction 


3.0 




3.0 




and rj veto 


1.0 


1.0 


1.0 


1.0 


Ks reconstruction 




0.7 




0.7 


Neural Net efficiency 


1.5 


1.0 


1.0 


1.0 


Fit Model 


0.7 


5.3 


2.9 


1.6 


Signal PDF bias 


0.9 


2.2 


1.6 


1.4 


Sum in quadrature 


3.5 


7.1 


5.3 


3.7 



5 RESULTS 

For the branching fraction calculation, we assume the production ratio, R^^^, is unity. R^^^ is 
defined as 

^+/o ^ r(r(4.) B+B-) 

T{T{As) ^ B^B^) ■ 

The measured branching fractions are shown in Table [2l The combined branching fractions are 
calculated from the sub- modes using the method of least squares, taking into account correlated 
systematic errors. 

To calculate the isospin asymmetry Aq-, we combine the branching fractions, the ratio of the 
B^ and B^ lifetime t+/to, and the production ratio R'^^^ according to 

Ao- = ^(/i?+/°^-l), (4) 

where I is 
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- B{B*- K*--f)' 

to obtain the isospin asymmetry 

Ao- = 0.029 ± 0.019 ± 0.016 ± 0.018. 

The first and second errors are statistical and systematic, respectively. The last error comes from the 
error on the production ratio, and we have used t+/to = 1.071 ±0.009 [21], = 1.020±.034 [22]. 
In addition, to obtain Eq. HI we have used the approximation that /, i?^/^, and t^/tq are all close 
to unity. The 90% confidence interval for Aq- including systematic uncertainties is 

-0.021 < Ao- < 0.079. 
The corresponding time- integrated CP asymmetry (table [2]) is 

A = -0.009 ±0.017 ±0.011, 
while the 90% confidence interval for A is 

-0.043 <A< 0.025. 

The combined asymmetries are calculated using the same method as the branching fractions. 

To ensure that we are measuring real K* mesons in data, we widen the K* mass selection to 
be 0.7 < rriK-K < 1-1 GeV/c^, refit the data, and make an sPlot [23] of the K* mass. We then fit a 
relativistic P-wave Breit-Wigner line shape to the sPlot. This is shown in Figure [5l We combine 
the measurements of mean of the K* mass and the width for the charged and neutral mesons 
separately to obtain the results in Table HI The results are consistent with the PDG values. 





Data 


PDG Value 


K* meson 


m (MeV) 


r(MeV) 


m (MeV) 


r(MeV) 




894.34± .63 


47.1± 1.4 


896.00± .25 


50.3± 0.6 


K*+ 


892.88± .80 


46. 7± 1.8 


891.66± .26 


50.8± 0.9 



Table 4: The combined results of the fits to the rriKiT spectrum shown in Figure [H Also shown are 
the PDG values. 



6 CONCLUSIONS 

We present a preliminary measurement of the branching fractions B{B^ K*^^) = (4.58 ± 0.10 ± 
0.16) X 10-5 and B{B+ K*+j) = (4.73 ± 0.15 ± 0.17) x 10"^. We use these results to calculate 
the isospin asymmetry at the 90% confidence interval to be —0.021 < Aq- < 0.079. We also present 
a preliminary measurement of the time-integrated CP asymmetry at the 90% confidence interval 
to be —0.043 < A < 0.025. These results are all improvements over previous measurements, as well 
as being consistent with SM expectations. 
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<m^^^> = 0.89447 ± 0.00068 

r= 0.0491 ±0.001 6 
^ X^/ndf = 1.177 
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Figure 5: Relativistic P-wave Breit-Wigner line 
of the sPlot of data for the a) K*^ K+tt', 
K*+ KsT^^ modes. 
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(d) K*-^ KsiT+ mode 

shape fit to the Kit invariant mass distribution 
b) K*° KsTT^, c) K*+ K+TT^, and d) 
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